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Rotational transition lines of CO play a major role in molecular radio astronomy and in particular in the study of star formation and the Galactic 
structure. Although a wealth of data exists in the Galactic plane and some well-known molecular clouds, there is no available CO high sensitivity 
• all-sky survey to date. Such all-sky surveys can be constructed using the Planck HPI data because the three lowest CO rotational transition lines 
CO at 115, 230 and 345 GHz significantly contribute to the signal of the 100, 217 and 353 GHz HPI channels respectively. Two diff'erent component 
separation methods are used to extract the CO maps from Planck HPI data. The maps obtained are then compared to one another and to existing 
CO external CO surveys. Prom these quality checks the best CO maps in terms of signal to noise and/or residual foreground contamination are selected. 
Three sets of velocity-integrated CO emission maps are produced: type 1 maps of the CO /=1^0, 7=2^1, and 7=3^2 rotational transitions with 
^ low foreground contamination but moderate signal-to-noise ratio; type 2 maps for the 7=1^0 and 7=2^1 transitions with a better signal-to-noise 
. ^ ratio; and one type 3 map, a line composite map with the best signal-to-noise ratio in order to locate the faintest molecular regions. The maps are 
described in detail. They are shown to be fully compatible with previous surveys of parts of the Galactic Plane and also of fainter regions out of the 
^ Galactic plane. The Planck HPI velocity-integrated CO maps for the 7= 1 ^0, 1 , and 7=3 -^2 rotational transitions provide an unprecedented 

all- sky CO view of the Galaxy. These maps are also of great interest to monitor potential CO contamination on CMB Planck studies. 
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ABSTRACT 



1. Introduction 

This paper, one of a set associated with the 2013 release of 
data from the Planck^ mission (Planck Collaboration I 2013), 



^ Planck (http://www.esa.int/Planck) is a project of the 
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describes the construction and validation of full- sky carbon 
monoxyde (CO) maps from Planck data. 

The interstellar medium (ISM) represents about 10 - 15% of 
the total mass of the Milky Way. The ISM is a mixture of atomic 
and molecular gas, the latter containing around 50% of its mass 
whilst filling only a tiny fraction of the volume (see Ferriere 
2001; Cox 2005, for a general introduction and review). The 
cold neutral gas is confined close the Galactic disk mid-plane 
(about 100 pc scale height), in clouds with varying molecular- 
to-atomic ratios. Molecular clouds, where hydrogen is molecu- 
lar, are the sites of star formation, and as such, play a pivotal role 
in the interstellar matter cycle. Giant molecular clouds (GMCs) 
are the largest self-gravitating structures in spiral galaxies such 
as the Milky Way, and are also the most massive entities, reach- 
ing several times 10^ Mq, placing them at the top of the mass 
spectrum. Molecular clouds were discovered via the rotational 
emission line /=1^0 of carbon monoxide in its fundamental 
electronic and vibrational levels (Wilson et al. 1970; Penzias 
et al. 1972). Contrary to the atomic component of the neutral 
ISM, which is directly observable via the spin-flip HI All cm 
line, the bulk of molecular hydrogen is not directly observable 
in molecular clouds. Because CO is abundant, easily excited by 
collisions with H2, and easily observable from the ground, it is 
considered a good tracer of the molecular component of the ISM. 
It is also a dominant coolant of molecular gas. The ability to de- 
tect the CO /=1^0 line from the ground allows large surveys 
to be performed (Dame et al. 2001), leading to the observational 
evidence that the molecular gas is structured in clouds that are 
turbulent and harbour the formation of all stars in the Galaxy. 
The astrophysical importance of carbon monoxide could hardly 
be over-appreciated. 

Large-scale surveys the /=1^0 line of ^^CO, but also of 
^^CO and C^^O isotopologues, have been carried out with meter- 
sized radio telescopes, mostly through the fundamental rota- 
tional transitions. The most complete CO(l^O) survey is that 
of Dame et al. (2001), which covers the Milky Way at Galactic 
latitudes \b\ < 30°, with an efl'ective spatial resolution of 0.5°. 
CO and isotopologues were observed with the 4-m NANTEN 
telescope providing spectral maps at slightly higher spatial reso- 
lution, towards specific GMCs (Mizuno & Fukui 2004). In addi- 
tion, there exist a wealth of smaller CO(l^O) line surveys, such 
as for example in Orion and Monoceros (Wilson et al. 2005). 
Observations of CO, for the Orion and Monoceros regions, 
by Magnani et al. (1985); Hartmann et al. (1998); Magnani 
et al. (2000) have revealed the existence of molecular clouds at 
Galactic latitudes up to 55°. However, these high-latitude obser- 
vations provide only a limited view of the \b\ > 30° sky. 

High-resolution full-sky surveys of higher / CO transi- 
tions have never been carried out, essentially because these are 
much more time consuming than /=1^0 observations (atmo- 
spheric transmission is poorer, amd higher- spatial resolution re- 
quires finer spatial sampling). High-/ lines are expected to probe 
molecular gas with stronger excitation conditions (high density 
and/or warmer) better than does the 1 ^0 line. Combined with 
CO(l^O), observations of / > 2 lines would therefore provide 
global constraints on the physical conditions in the molecular 
ISM. Only specific regions have been mapped in the /=2^1 
CO line, e.g., the Galactic Centre region has been observed by 
Sawada et al. (2001). Maps of W3 and W5 (Bieging et al. 2010; 
Bieging & Peters 2011) are also available. For the /=3^2 line. 
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Fig. 1. The average spectral response for each of the HEI fre- 
quency bands. The vertical bars represent the CO rotational tran- 
sitions. 



the Galactic centre has been observed by Oka et al. (2012), Orion 
by Ikeda et al. (1999) and W3 by Bieging et al. (2010). 

The ESA Planck satellite^ was launched on the 14 May 
2009 in order to (primarily) measure with unprecedented pre- 
cision the temperature and polarization anisotropics of the cos- 
mological microwave background (CMB). It observes the sky in 
nine frequency bands covering 30-857 GHz with high sensitiv- 
ity and angular resolution from 3r to 5'. The Low Frequency 
Instrument (LEI; Mandolesi et al. 2010; BersanelH et al. 2010; 
MenneUa et al. 2011) covers the 30, 44, and 70 GHz bands with 
amplifiers cooled to 20 K. The High Frequency Instrument (HEI; 
Lamarre et al. 2010; Planck HEI Core Team et al. 2011a) covers 
the 100, 143, 217, 353, 545, and 857 GHz bands with bolome- 
ters cooled down to 0.1 K. Planck'^ sensitivity, angular reso- 
lution, and frequency coverage make it a powerful instrument 
for Galactic and extragalactic astrophysics as well as cosmology 
(Planck CoUaboration I 2013). 

The first seven CO rotational transition lines lie within the 
spectral bands of the HEI instrument. Of these, the first three, 
/=1^0, /=2^1, and /=3^2 at 115, 230 and 345 GHz, re- 
spectively, present the largest transmission coefficients making 
them a significant foreground component in the Planck intensity 
maps. In this paper, we extract full-sky CO maps for these three 
lines from the LEI and HEI data using component separation 
methods. The Planck intensity maps are presented in Sect. 2. In 
Sect. 3, we provide a brief description of HEI bandpasses and 
estimate the CO transmission coefficients for the most impor- 
tant rotational lines. Specifically tailored component separation 
methods for CO extraction in the Planck maps are detailed in 
Sect. 4. The Planck CO maps obtained using the above meth- 
ods are presented in Sect. 5. Uncertainties and foreground con- 
tamination on those maps are discussed in Sect. 6. The inter- 
nal validation of the Planck CO maps is presented in Sect. 7. 
Detailed comparison to existing external CO surveys is pre- 
sented in Sect. 8. Finally, we discuss the results and draw con- 
clusions in Sect. 9. 



^ Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions from NASA (USA) and 
telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 
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2. Planck data 

This paper is based on the Planck first 15.5 month survey mis- 
sion (two full-sky surveys) and uses the full-sky maps of the nine 
Planck frequency bands, and also the 100, 217 and 353 GHz 
full-sky bolometer maps. These maps are provided in HEALPix 
pixeHzation (Gorski et al. 2005) with A^side = 2048 at full res- 
olution. We refer to (Planck Collaboration II (2013); Planck 
Collaboration V (2013); Planck Collaboration VI (2013); ?); 
Planck HFI Core Team et al. (2011b); Zacchei et al. (2011) for 
the generic scheme of TOI processing and map-making, as well 
as for calibration. These Planck maps are given in /^cmb units, 
i.e. in temperature units referred to the CMB blackbody spec- 
trum. The scanning strategy of Planck consists of circles on the 
sky that correspond to diff'erent positions of the satellite spin 
axis. The latter is changed by 2.5' every 40 to 60 minutes in 
order to cover the full sky in about seven months. The data set 
acquired for each position of the satellite spin axis is called a ring 
and consist of 40 to 60 observations of the same circle on the sky. 
Thus, a noise map can be obtained for each frequency band or 
bolometer map from the diff'erence of the maps of the first and 
second half of the rings. The resulting noise maps are basically 
free from astrophysical emission and thus a good representation 
of the statistical instrumental noise and systematic errors. In the 
following we assume that the beam pattern of the maps can be 
well represented by eff'ective circular Gaussians with FWHM of 
32.^24, 27.^01, 13.^25, 9.'65, 7.'25, 4.^99, 4.^82, 4.^68, and 4.33 at 
30, 44, 70, 100, 143, 217, 353, 545, and 857 GHz, respectively 
(Planck Collaboration IV 2013; Planck Collaboration VII 2013). 



3. CO contribution in the Planck HFI channels 

3.1. HFI Spectral Response and CO Emission 

In order to isolate the narrow CO features from the remainder 
of the components (e.g., CMB, dust, etc.), precise knowledge 
of the instrument spectral response as a function of frequency 
is required. The original spectral resolution requirement for the 
spectral response of a given HFI detector was about 3 GHz; this 
corresponds to a velocity resolution of around 8000 km/s for 
the CO(l^O) line. As HFI does not have the ability to mea- 
sure spectral response within a frequency band during flight, the 
ground-based Fourier Transform Spectrometer (FTS) measure- 
ments provide the authoritative data on the HFI spectral trans- 
mission. With good S/N, spectral information can be inferred to a 
fraction (i.e., ~l/10th) of a spectral resolution element (Spencer 
et al. 2010). Thus, the pre-flight calibration FTS measurements 
(Pajot et al. 2010), which were carried out at a spectral resolution 
of about 0.6 GHz, may be used to estimate the spectral response 
at a resolution equivalent to 50 km s"^ for the CO(l^O) line. 
For spectral regions near CO rotational transitions, therefore, the 
spectral response was oversampled by a factor of around 10 us- 
ing an interpolation based on the instrument line shape of the 
FTS. Table 1 lists the relevant CO transitions, and the frequency 
ranges over which the HFI spectral response was oversampled 
(Planck Collaboration IX 2013). The oversampling ranges were 
extended to include all of the common CO isotopologues. The 
vertical bars above the spectral response curves of Fig. 1 illus- 
trate these oversampled regions. 

3.2. Spectral Band CO conversion coefficients 

Unit conversion coeflftcients between CO brightness tempera- 
ture and CMB temperature are determined in much the same 



CO 




2 ' ' ' 4 ' ' ' 6 ' ' ' 8 
Bolometer pair number 



CO 




00 ^ 



2 4 6 8 

Bolometer pair number 

Fig. 2. Comparison of normalised CO conversion coeflScients 
computed from the bandpass measurements (black) or estimated 
on the sky (blue) for the 100 (top), 217 (middle) and 353 GHz 
(bottom) channels. See Table 2. 

manner as the other HFI unit conversion coefl&cients (see Planck 
Collaboration IX 2013): 

'''''''' J dv^I^^^ ^ ^ 

where is the spectral transmission of bolometer b at fre- 
quency y, and 1^^ and 7™^ are the CO and CMB intensities. 
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Table 1. Rotational ^^CO transitions within the HFI bands and 
over- sampled regions. 



Band CO transition vq Oversampling 

[GHz] [/upper ^ /lower] [GHz] [GHz] 

100 1^ 115.271 109.6- 115.4 

217 2^ 1 230.538 219.3- 230.8 

353 3 ^ 2 345.796 329.0- 346.2 

545 4 ^ 3 461.041 438.6- 461.5 

545 5 ^ 4 576.268 548.3- 576.8 

857 6 ^ 5 691.473 657.9- 692.2 

857 7 ^ 6 806.652 767.5- 807.5 

857 8 ^ 7 921.799 877.0- 922.7 

857 9 ^ 8 1036.912 986.6- 1037.9 



respectively. The CO velocity-integrated emission (VIE) is ex- 
pressed as the product of Ray leigh- Jeans temperature and spec- 
tral line width in velocity units, i.e., Krj km s"^ For the lowest 
rotational CO transitions a Doppler shifted line profile may be 
assumed with v = vcoCl + v/c)"^ ^ vco(l - v/c) for v «c c. As 
the CO transitions occur at discrete frequencies, with a Doppler 
line- width much less than the transition frequency (about 10^ Hz 
cf. 10^^ Hz), and much narrower than the available knowledge of 
the HFI detector spectral response (~ 10^ Hz), the velocity inte- 
gration may be approximated by a delta function distribution at 
vco, i.e. ^vco- 

Bandpass CO conversion coefficients for the first three CO 
transitions were determined using the above relation (see Planck 
Collaboration IX 2013)) ; similar data are available for the in- 
dividual HFI detectors and for the other transitions within HFI 
bands. Table 2 (first column) lists the relative values of these 
CO conversion coefficients for the ^^CO isotopologue (aver- 
aged for each a and b pair of polarized bolometers) when nor- 
malized to the average transmission of all bolometers in the 
band (this format is chosen to ease comparison to the sky- 
calibrated conversion coeflftcients, see next section). To recover 
the physical conversion factor, values in the table must be mul- 
tiplied by the average of ^^CO conversion coefl&cients, namely, 
1.42 X IQ-^KcMB km s'^ at 100 GHz, 4.50 x IQ-^Kcmb km s'^ 
at 217 GHz and 17.37 x IQ-^Kcmb km s'^ at 353 GHz. 

As will be described in Sect. 6.2.2, the ^^CO isotopologue 
may contribute to the Planck CO maps. The conversion coef- 
ficients of ^^CO are computed in the same way as described 
for i^co, but for the ^^CO transitions at v^^'^^ = 110.2 GHz, 

y[f"^^ = 220.40 GHz, and v^^'^^ = 330.6 GHz. The average con- 
version factor over all bolometers in the 100, 217, and 353 GHz 
channel are 1.62 x 10"^Kcmb km s"\ 3.63 x 10"^Kcmb km s"\ 
and 13.0 X 10"^Kcmb km s"\ respectively. 

3.3. Sky-calibrated CO conversion coefficients 

In the following sections we will use CO conversion coefficients 
to extract maps of the CO emission from the Planck frequency 
maps. For this purpose, we need an accurate estimation of the 
relative CO conversion coefficients between bolometers (see ap- 
pendix B for further discussions). The bandpass CO conversion 
coefficients present two main problems. First, the estimation of 
the spectral band transmissions of each of the bolometers may be 
aff'ected by systematic errors as we are dealing with narrow lines. 
Indeed, HFI bandpasses are not sampled to a sufficient resolution 



to allow for a satisfactory CO extraction using these coefl&cients 
and then need to be interpolated as discussed above. Second, 
both ^^CO and ^^CO (and other isotopologues) are present in the 
Planck maps. The emissions from ^^CO and ^^CO are spatially 
correlated with a varying ratio across the sky and we can not dis- 
criminate between the two. Thus, an accurate determination of 
the relative CO conversion coefficients between bolometers of 
the same Planck channel is difficult when using bandpass infor- 
mation only. 

A way around these issues has been found in estimating the 
CO conversion coefficients directly from sky measurements of 
well-known molecular clouds. For these regions, we can assume 
- to first order - that the map of the sky emission for a detec- 
tor is a linear combination of CO and thermal dust emissions 
weighted by the CO and thermal dust conversion coefl&cients of 
the detector. Thus, using an external CO emission template, Pqo, 
obtained from the Dame et al. (2001) ^^CO /=1^0 survey and 
a dust emission template, the 545 GHz Planck-HFl channel map 
(in KcMB units), it is possible to determine sky-calibrated CO 
and dust transmission coefficients F^^^ ^ and F^^^^ ^ of bolometer 
b. Notice that we implicitly assume here a perfect spatial corre- 
lation between the ^^CO and ^^CO components. 

Neglecting other astrophysical components (i.e., QMB, free- 
free, synchrotron, AME and point sources), we perform a simple 
linear fit between each bolometer map in the native units 
(Kcmb) and the CO and dust templates: 

= F^o,b ^CO + ^dmt,/. , (2) 

where F^^ ^ and F^^^^ ^ are the results of the linear regression. To 
avoid contamination from dust polarization emission we com- 
bine the maps of pairs of polarized bolometers into a single map 
and a single CO conversion coefficient is computed for the pair 
of detectors. 

Due to the spatial correlation between the CO and dust emis- 
sion, we expect the CO conversion coefficients, ^^0/?' 
fected by dust contamination and so biased. We simply assume 
here that the estimated CO conversion factor can be expressed as 
^co b ~ ^ ^co b^P corresponding to an offset and a shift in total 
amplitude. However, we will consider that thermal dust conver- 
sion coefficients are not aff'ected by CO contamination signifi- 
cantly and thus i^dust,/. = ^L,/.- 

To compute a and p we solve the following system of equa- 
tions. First, we construct dust weighted bolometer (pairs of 
bolometers) maps. Then, we obtain a first set equations by im- 
posing that the difference of two dust weighted bolometer (pairs 
of bolometer) maps, 

K Ml, 

rfit p rfit p 

must be correlated with the CO template (weighted by the dust 
map). A second set of equations is obtained by searching for the 
factor y that minimizes the correlation of the difference 

Ml Ml 

^ rx ^— 

rfit p ' rfit p 

with the CO dust weighted template, Pco/Pd- Uncertainties in 
the final CO conversion coefl&cients are obtained from the dis- 
persion of the coefficients found for different sky regions. For 
the CO(l^O) line we use sky regions for which the Dame et al. 
(2001) map is above 2 Krj km s"^ For the /=2^1 and /=3^2 
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Table 2. Relative CO conversion coefficients of HFI bolometers 
(pairs of bolometers) in the 100, 217 and 353 GHz channels nor- 
malized to the average conversion coefficient for each channel. 



Bolo ID 


77BP 


/7BP 




/=1^0 


100-1 (a+b)/2 
100-2 (a+b)/2 
100-3 (a+b)/2 
100-4 (a+b)/2 


0.82 + 0.10 
0.94 + 0.09 
0.99 + 0.11 
1.24 + 0.10 


1.03 + 0.12 
0.97 + 0.10 
0.87 + 0.14 
1.13 + 0.24 


0.83 + 0.01 
1.09 + 0.01 
0.88 + 0.01 
1.19 + 0.01 


/=2^1 


217-1 
217-2 
217-3 
217-4 

217-5 (a+b)/2 
217-6 (a-Hb)/2 
217-7 (a+b)/2 
217-8 (a+b)/2 


1.10 + 0.03 
0.94 + 0.05 
1.13 + 0.12 
1.05 + 0.06 
0.97 + 0.05 
0.88 + 0.06 
0.97 + 0.15 
0.96 + 0.05 


0.98 + 0.06 
0.94 + 0.13 

1.07 + 0.06 
0.88 + 0.20 
0.89 + 0.06 
1.06 + 0.12 

1.08 + 0.05 
1.08 + 0.04 


1.03 + 0.01 
0.95 + 0.01 
1.16 + 0.01 
1.14 + 0.01 

1.00 + 0.01 
0.79 + 0.01 

1.01 +0.01 
0.92 + 0.01 


7=3^2 


353-1 
353-2 

353-3 (a+b)/2 
353-4 (a+b)/2 
353-5 (a+b)/2 
353-6 (a+b)/2 
353-7 
353-8 


0.98 + 0.04 
1.00 + 0.05 
1.11 +0.07 
0.90 + 0.09 
0.89 + 0.06 
0.91 +0.14 
1.14 + 0.07 
1.06 + 0.09 


0.66 + 0.07 
1.05 + 0.05 
1.22 + 0.11 
1.13 + 0.08 
1.24 + 0.09 
0.89 + 0.11 
0.12 + 0.04 
0.80 + 0.08 


0.94 + 0.01 
0.95 + 0.01 
1.07 + 0.01 
0.84 + 0.01 
0.96 + 0.01 
1.01 +0.01 
1.15 + 0.01 
1.07 + 0.01 



survey that it is used as a CO template to the original ^cmb units 
of the Planck HFI maps. To be useful to scientific purpose these 
maps must be converted into the emission of the transition line. 
This is done by recalibrating the Planck CO maps so that they 
are in units of the expected ^^CO contribution at the transition 
frequency: 



y^/final _ y^sky 



Here M'^^ is the CO map of any transition calibrated on the 
Dame et al. (2001) data, using the sky CO conversion coeffi- 
cients. The quantities {F^^) and {F^qq) are the average across 
detectors of the bandpass and sky ^^CO conversion coefficients. 



(3) 



4. Extracting CO from Planck data 

To extract the CO emission from the Planck maps, three main 
approaches have been considered: i) a single channel analysis 
(Sect. 4.2.1), ii) a multi-channel approach (Sect. 4.2.2) and fi- 
nally iii) a multi-line approach (Sect. 4.2.3) using fixed CO line 
ratios. While the two first methods allow us to reconstruct spe- 
cific CO transition lines, the third one does not discriminate be- 
tween CO transitions, but yields the best signal-to-noise ratio. 
Before detailing the specifics of these approaches, we present in 
Sect. 4.1 the various component separation methods that have 
been used and adapted to the specifics of CO observations with 
the Planck satellite. 



lines we consider only the nine brightest CO clouds in the Dame 
et al. (2001) map. 

The sky CO conversion coefficients^ (normalized to their 
mean value) obtained from this analysis are presented in the 
right column of Table 2 and compared to ^^CO bandpass con- 
version coefficients in Fig. 2. The discrepancy between the rel- 
ative values of the bandpass and sky determined CO conversion 
can be explained from the quoted uncertainties. A more detailed 
comparison of sky and bandpass determined CO conversion co- 
efficients, and a description of the systematic uncertainties is 
given in the Planck Collaboration IX (2013) companion paper. 
However, it is important to notice that they are not expected to 
be equal as the sky-based CO conversion coefficients account 
also for the contribution of other spatially correlated CO iso- 
topologues, mainly ^^CO and other molecular lines. This is ex- 
tensively discussed in Sect. 6.2.2. 

Finally, note that for the /=2^1 and /=3^2 lines the maps 
obtained from the sky-calibrated CO conversion coefficients are 
calibrated in the Dame et al. CO(l^O) line units and need to be 
re-calibrated to their actual frequencies; this is the purpose of the 
next section. 

3.4. Unit convention 

In the following the Planck CO maps are extracted using the 
sky determined CO conversion coefficients. As discussed above 
these coefficients convert from the units of the Dame et al. (2001) 

^ Compatible results within statistical errors are obtained when per- 
forming a similar, but independent, analysis with Commander. This re- 
sults shows these conversion coefficients are robust against some possi- 
ble systematic effects on the method. 



4.1. Component separation methods 

Component separation algorithms as presented in Planck 
Collaboration XII (2013) are mainly specifically tailored for 
CMB extraction and its statistical analysis. Thus, for this paper 
we have adapted and tested several algorithms for CO extraction 
from which we have selected two. Here we give a brief overview 
of them and of their main characteristics, while a more detailed 
description can be found in the references given below. 

4.1.1. MILCA 

The MILCA (Modified Internal Linear Combination Algorithm) 
method (Hurler 2012) was specifically developed within the 
Planck collaboration or the reconstruction of thermal Sunyaev- 
Zeldovich effect and CO contributions on the Planck maps. 
MILCA is an extension of the standard ILC algorithm origi- 
nally aiming at CMB extraction (Bennett et al. 2003; Eriksen 
et al. 2004). MILCA provides a flexible way of selecting wanted 
and un- wanted spectral components (Hurler et al. 2010; Hurler 
2012) and corrects for the noise bias in standard ILC algo- 
rithms. Furthermore, the weights of the internal linear combina- 
tion can be computed both in real and harmonic space to improve 
the component separation efficiency. MILCA was tested against 
CO-oriented simulations, showing no bias in the reconstruction 
when the CO conversion coefficients are perfectly known (see 
Appendix B). 

4.1.2. Ruler 

Ruler is an inversion algorithm, which provides the general- 
ized least squares solution (GLSS), given a parametric model of 
the Galactic emission. In a nutshell, input data, either bolometer 
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Type 1 J=1-0 CO map Type 1 J=1-0 stddev 





Type 3 J=1-0 CO map 



Type 3 J=1-0 stddev 




-5.0^^^^^^^^^^^^ .^^^^^^^^^^^H 30.0 0.37^^^^^^^^^ .^^^^^^^^H 3.5 

K Km/s K Km/s 

Fig. 3. Full-sky CO /=1^0 maps (left) and their respective statistical error maps (right) for the three types of Planck CO products. 
For display purposes only, the type 1 and type 3 maps are smoothed to 15' to ease the comparison to the type 2 maps. The error 
maps correspond to the non-smoothed product. The maps are in Galactic coordinates and follow the HEALPix pixelization scheme 
forA^side = 2048. 



maps or frequency maps, are linearly combined, on a pixel-by- 
pixel basis, according to weights that account for both the fore- 
ground spectral properties and the Planck channel specifications, 
namely instrumental noise and bandpass. A comprehensive de- 
scription of the method is given in Planck Collaboration XII 
(2013). Validation on simulations are presented in Appendix B. 



Note that Ruler and MILCA have both been tested in the 
single- and multi-channel configurations described in Sect. 4.2 
and were found to produce compatible results. This brings con- 
fidence in the robustness of these independent algorithms. 



4.1.3. Coimnander -Ruler 

The Coimnander-Ruler component separation pipeline consists 
of two steps: i) the Planck channel maps are brought to a com- 
mon resolution and the likelihoods of the non linear degrees of 
freedom of a chosen parametric foreground model are jointly 
sampled, through a Gibbs sampling Monte-Carlo algorithm 
(Coimnander, Eriksen et al. 2008); ii) for each sample of the de- 
rived distribution, the linear degrees of freedom, namely the am- 
plitude maps of the components, are computed at the full Planck 
resolution via a generalized least squared solution (Ruler, de- 
scribed above). The first and second moments of the posterior 
distribution of a parameter define its mean value and uncertainty, 
which accounts for both the instrumental and foreground mod- 
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Type 1 J=2-1 CO map Type 1 J=2-1 stddev 





Fig. 4. Full-sky CO J=2^1 maps (left) and their respective statistical error maps (right) for type 1 (top) and type 2 (bottom) Planck 
CO products. The type 1 map has been smoothed to 15' for display purposes but the error map corresponds to the non smoothed 
product. The maps are in Galactic coordinates and follow the HEALPix pixelization scheme for A^side = 2048. 



Type 1 J=3-2 CO map 



Type 1 J=3-2 stddev 




Fig. 5. TYPE 1 full-sky CO J=3^2 maps (left) and its statistical error map (right). The map has been smoothed to 15' for display 
purpose but the error map corresponds to the non smoothed product. The maps are in Galactic coordinates and follow the HEALPix 
pixelization scheme for A^side = 2048. 



elling error. Coimnander- Ruler has been extensively tested us- 
ing the FFP6 simulations (see Planck Collaboration ES 2013) as 
part of Planck component separation, a thorough description of 
which is provided in (Planck Collaboration XII 2013). 



4.2. CO extraction strategies 

4.2.1 . Single-channel approach - type 1 maps 

In the single-channel approach we exploit differences in the 
spectral transmission of a given CO line among the bolometers 
(pairs of bolometers) of the same frequency channel. The main 
advantages of this solution are to give access to the first three 
transitions of CO at the native resolution of the Planck maps and 
to avoid contamination from other channels. However, this type 
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Table 3. Main characteristics of the Planck CO products. 



Type 




Line 


Resolution 

n 


Noise / 15'pixel 
[Kkms-i] 


side 


Method 


Components considered 


Data used 
Freq [GHz] 


TYPE 


1 


CO(l^O) 


9.65 


1.77 


2048 


MILCA 


CO, CMB 


100 (bolo. maps) 


TYPE 


1 


C0(2^1) 


4.99 


0.74 


2048 


MILCA 


CO, CMB, dust 


217 (bolo. maps) 


TYPE 


1 


CO(3^2) 


4.82 


0.73 


2048 


MILCA 


CO, dust 


353 (bolo. maps) 


TYPE 


2 


CO(l^O) 


15 


0.45 


2048 


Ruler 


CO, CMB, dust, free-free 


70, 100,143,353 


TYPE 


2 


C0(2^1) 


15 


0.12 


2048 


Ruler 


CO, CMB, dust, free-free 


70,143,217,353 


TYPE 


3 




5.5 


0.16 


2048 


Commander -Ruler 


CO, CMB, dust. 


30-857 



power-law sync/free-free 



of CO extraction results in a lower signal-to-noise ratio due to 
the use of individual bolometer (pair of bolometers) maps. 

The CO /=1^0, J =2^1 and /=3^2 maps are obtained 
with MILCA, using all bolometers in the 100, 217 and 353 GHz 
channel, respectively, and are denoted the Planck type 1 CO 
product. Different constraints are applied depending on the line 
under scrutiny: 

- The /=1^0 map is obtained by requesting cancellation of 
a flat spectrum (i.e. the CMB) while preserving CO (using 
the sky-calibrated conversion coefficients given in Table 1). 
Notice that in this case the diff'use Galactic foreground con- 
tamination can be well approximated by a flat spectrum 
within the band and there is no need for extra constraints. 

- At 217 GHz, dust becomes more of a major foreground and 
should be dealt with. The /=2^1 map is therefore extracted 
using an extra constraint on the dust transmission in the dif- 
ferent bolometers, computed from a grey-body spectrum de- 
fined as oc v^By(Td) where T = UK, (3=1.6 and B^iT) 
is the Planck function. The choice of this dust spectrum is 
discussed in Sect. 6.2.1. 

- At 353 GHz, the CMB becomes sub-dominant compared to 
the dust and we only request keeping the CO and removing 
the dust. We also found that the dust spectrum used for the 
/=2^1 line is not optimal at 353 GHz and thus we used 
the sky-based dust conversion coefficients F^^^^ which were 
also fitted in the linear regression of Eq.(2) and for which we 
obtained better results. 

4.2.2. Multi-channel approach - type 2 maps 

The multi-channel approach makes use of the intensity maps in 
several channels to isolate the CO contribution from the trans- 
mission of CO, CMB, dust and free-free emission in the diff'er- 
ent channels. The Ruler method was selected to be used in this 
configuration, to extract the /=1^0 and /=2^1 lines. 

In order to construct the CO(l^O) line, we use LFI 70 GHz, 
HFI-100, 143 and 353 GHz channels, while C0(2^1) is ob- 
tained using LFI 70 GHz, 143, 217 and 353 GHz channels. AH 
maps are smoothed to 15'. To solve for the two CO transitions, 
we assume that the CO contribution to the 353 GHz channel is 
negligible"^. The CO(l^O) and C0(2^1) maps obtained using 
the multi-channel method constitute the Planck type 2 CO prod- 
uct. 

For both lines, the requirement is to extract CO while fitting 
at the same time for CMB, dust and free-free emission. By con- 
struction, the CMB spectrum is flat across channels. For the dust, 

^ While this assumption eases the CO extraction, it also leads to an 
overall systematic calibration error that can be corrected for afterwards 
(see Sect. 6.2). 



the same grey-body spectrum as for the single-channel /=2^ 1 
line is assumed. Free-free emission is modelled as a power-law 
spectrum oc y"^-^^. The dust and free-free transmission coeflft- 
cients are obtained by integrating their spectra over Planck'^ 
bandpasses. 

4.2.3. Multi-line approach - type 3 map 

While the multi-channel approach provides single line maps 
at relatively high signal-to-noise ratio, it is possible to use a 
multi-line approach to increase even further the S/N ratio, po- 
tentially allowing the discovery of new faint molecular clouds 
at high Galactic latitudes. For this, we need to assume that the 
line ratios (CO(2^1)/CO(1^0) and CO(3^2)/CO(1^0)) are 
constant across the sky. This approach is used as part of the 
Coimnander-Ruler pipeline employed for the component sep- 
aration. 

The TYPE 3 CO map results from a seven-band run of this 
pipeline, including 30 to 353 GHz channel maps, with the sky 
modelled as a superposition of CMB, CO, dust - treated as 
a modified blackbldy - and a power law to describe the low- 
frequency Galactic emission. The dust optical depth and the dust 
temperature, as well as the low-frequency component spectral 
index are fitted at every pixel. Since the total number of pa- 
rameters would exceed the number of frequencies considered, 
the three CO line emissions are assumed to be perfectly cor- 
related. A single CO map is solved for, the so-called type 3 
map, using the average CO bandpass transmission (see Sect. 3.2) 
in each channel, whereas the CO line ratios are given by the 
posterior average of the distribution obtained from a dedicated 
Coimnander run on small bright CO regions (Taurus, Orion, 
Polaris, etc.). The line ratios found for CO(2^1)/CO(1^0) and 
CO(3^2)/CO(1^0) are 0.595 and 0.297, respectively. Notice 
that there are significant spatial variations on these ratios as dis- 
cussed in Sect. 6.2.2. For the analysis we have selected bright 
CO regions with large ratios. 



5. Planck CO maps 

The Planck CO delivery consists of three types of products, cor- 
responding to the three extraction methods described above: 

- TYPE 1 maps are extracted using the single-channel approach 
described in Sect. 4.2.1 and come at the native resolution of 
the corresponding Planck-Y{F\ channel; 

- TYPE 2 maps come from the multi-channel method given in 
Sect. 4.2.2 and have a resolution of 15'.; 

- The TYPE 3 map comes from the multi-line approach de- 
scribed in Sect. 4.2.3, and note that it assumes fixed CO line 
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ratios in order to obtain the highest possible signal-to-noise 
ratio, with is delivery at a resolution of ~ 5'.5. 

The main characteristics of all Planck CO maps are gathered in 
in Table 3 while details of the Planck CO product is given in 
Appendix A. The type 1, type 2 and type 3 terminology will be 
used throughout the paper. All maps discussed below are in units 
of Krj km s"^ at the transition line they represent. 

5.1. The CO J=1^0 line at 1 15 GHz 

The TYPE 1 (top), TYPE 2 (middle) and type 3 (bottom) Planck 
CO(l^O) maps are given in the left column of Fig. 3. The right 
column corresponds their statistical error maps and will be dis- 
cussed further in Sect. 6. 

The resolution of the type 1 product is 9.65', i.e. the native 
resolution of the Planck-HFI 100 GHz channel. The type 2 prod- 
uct has a 15' beam due to additional smoothing required to com- 
bine several Planck channels: in particular the use of Planck- 
LFI 70 GHz channel is essential to remove the free-free emis- 
sion (see Sect. 6.2). Finally, the type 3 has a varying resolution 
across the sky, resulting from the Ruler solution, and a beam 
profile is computed from FFP6 simulations. A good Gaussian 
approximation can be achieved with a FWHM of 5.5'. A simple 
eye inspection finds a good overall agreement between the maps; 
an in-depth comparison is conducted in Sect.7. 

5.2. The CO J =2^ 1 line at 230 GHz 

The C0(2^1) line can only be extracted using the single- 
channel or the multi-channel approach. Therefore the 230 GHz 
Planck CO product consists of type 1 (5' resolution) and type 2 
(15' resolution) maps only, which are displayed in the left col- 
umn of Fig. 4. The dust emission increases with frequency and 
becomes more of an issue for the /=2^1 line extraction. This 
will be discussed at length in Sect. 6.2, but this issue is already 
clearly visible when comparing the two products: the type 2 CO 
presents a diff'use dust emission throughout the Galactic plane 
that is not present in the type 1 map. 



noise map by the square root of the number of hits in each pixel, 
namely 



1 



CTb = 



stddev 



[Ml -Ml) VA^ 



(4) 



where is the hit number map and Ml and are the first and 
last half-ring bolometer maps respectively. This standard devia- 
tion of each bolometer map CTb is then propagated quadratically 
to the CO map level using the weights found by MILCA for the 
linear combination. 

This cannot be achieved in such a straightforward fashion 
for the TYPE 2 maps that have been smoothed to 15', resulting in 
correlated noise. Given the size of the maps, a fully dense ma- 
trix description of the noise is not feasible. However, the set of 
1 ,000 realistic noise simulations has been processed through the 
pipeline (first smoothed and then linearly combined), and used 
to compute a cr-map, which is a good pixel noise approxima- 
tion. The TYPE 3 standard deviation map is obtained in a similar 
fashion. 

The (T-maps of all Planck CO maps are plotted in the right 
columns of Figs. 3, 4 and 5. We remind here that for the type 1 
and TYPE 3 CO(l^O) maps, the cr-maps correspond to the 
maps at their native resolutions, while the CO maps have been 
smoothed to 15' for display purpose. These standard deviation 
maps have all been validated by checking that their mean at high 
latitude was in agreement with standard deviation measured di- 
rectly in the CO map at these locations (where no signal is ex- 
pected, see lower panel of Fig. 17). 

Using a common resolution of 15' to compare the high- 
latitude noise level in the maps, we find the standard deviation of 
the CO(l^O) maps to be typically 1.77 Krj km s"^ for the type 1 
map, 0.45 Krj km s"^ for the type 2 map and 0.16 Krj km s"^ 
for the TYPE 3 map. These uncertainties can be compared with 
typical uncertainties on ground-based surveys for CO(l^O) for 
example 0.6 Krj km s"^ for the and 1.2 Krj km s"^ for the Dame 
et al. (2001) and NANTEN (Mizuno & Fukui 2004) surveys. At 
the same resolution, the type 1 and type 2 C0(2^1) maps have 
standard deviations of 0.74 Krj km s"^ and 0.12 KRJkms"^ re- 
spectively, while it is 0.73 Krj km s"^ type 1 CO(3^2) map. 



5.3. The CO J =3^2 line at 345 GHz 

As mentioned previously, the single-channel approach is the 
only way to extract this higher / CO line and the correspond- 
ing TYPE 1 CO map is shown in Fig. 5. This map has a resolution 
of4'.82. 



6. Uncertainties and foregrounds 

6.1. Statistical errors 

Statistical uncertainties in the CO maps may be obtained using 
so-called half-ring diff'erences^ (Planck HFI Core Team et al. 
2011b). For the type 1 product, standard deviation ((T-)maps 
are generated at the level of individual bolometer maps, starting 
from the half-ring diff'erences noise map that is whitened using 
the number of hits in each pixels. A given bolometer cr-map, 
CTb, is then obtained by dividing the standard deviation of the 



^ First and last ring sets are independent data sets built using the first 
and second half of the stable pointing periods (see Sect. 2). 



6.2. Foregrounds 

All Planck CO products suff'er from systematic eff'ects and 
foreground contamination that need to be characterized; point 
sources or emission from other CO lines aff'ect all Planck 
CO products while dust, free-free emission and the Sunyaev 
Zeldovich eff'ect (SZ) are important for type 2 and type 3 maps 
only (the type 1 product is the most immune to contamination 
as it relies on single channel information). Some of these fore- 
ground emissions (e.g., CMB, dust, free-free), need to be dealt 
with at the stage of the component separation. Others, like point 
sources and the thermal SZ eff'ect in clusters of galaxies, may 
simply be masked afterwards. We discuss each of these fore- 
grounds in detail below apart from CMB emission for which the 
electromagnetic spectrum is well-known and can be explicitly 
nullified within calibration errors in the Planck frequency maps, 
as discussed in Sect. 4.1. 

6.2.1. Dust 

Dust emission is the main foreground as far as CO is concerned 
and is often strong in the same regions, e.g. star forming regions. 
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Taurus 



: • TYPE 1 
: •TYPE 2 
I •TYPE 3 




10 20 30 40 

Dust [HE I 545 GHz(IVUy sf-')] 

Polaris 




Dust pHFI 545GHz(MJv sr^)] 



Fig. 6. Comparison of 3 types of CO product with a dust template 
map for the Taurus (top), Orion (middle) and Polaris (bottom) 
molecular clouds. 



The multi-channel approach (i.e. type 2 products) is particularly 
sensitive to the choice of the dust spectrum. One assumption we 
make in these maps is that the dust spectrum is constant over 
the sky and described by a grey-body with Tdust = 17 K and 
y^dust = 1-6. This assumption is not correct since both the spectral 
index and temperature of the dust are known to vary across the 
sky (see for example Planck Collaboration XIX 2011; Fauvet 
et al. 2013). These values are, however, a good representation of 



the dust found in CO-rich regions as can be seen on the y^dust and 
Tdust maps of Planck Collaboration XII (2013). 

First, let us remark that as far as the dust transmission in HFI 
is concerned, /5dust and Tdust are degenerate quantities. Testing 
different values of the /^dust index^, we find this dust model 
to be the best compromise for the multi-channel CO extrac- 
tion. At Tdust = 17 K, choosing y^dust = 1-5 results in too 
much dust removal and a negative residual in the Galactic plane. 
Conversely, y^dust = 1.7 generates a much larger remaining dust 
emission. In the outer regions of the type 2 maps (/ > 90° and 
/ < 270°), changing y^dust by ±0.1 yields an increase/decrease of 
< IKrj km s"^ at 115 GHz and < 2 Krj km s"^ at 230 GHz. 
Still, these are larger than statistical errors. The inner regions of 
the Galactic plane are far more sensitive to a change of the dust 
spectrum, with a shift of about 10 K km s"^ at 115 GHz and 
15 K km s"^ at 230 GHz, but the signal is much larger. 

Given our chosen dust model. Fig. 6 shows the correlation 
between the types of /=1^0 CO maps discussed above and 
Planck'^ 545 GHz map used here as dust template. All maps 
were smoothed to 30'. The comparison has been done in three 
molecular clouds hosting diff'erent environments: i) Orion, a very 
active massive star forming region; ii) Taurus, that hosts low- 
mass star formation (Kenyon et al. 2008); and iii) Polaris, a 
high-Galactic latitude translucent cloud with little to no star for- 
mation and presenting both atomic and molecular gas (Miville- 
Deschenes et al. 2010). Figure 7 presents smoothed images of 
the three regions for the type 1 and type 2 Planck CO maps 
and for Planck'^ 545 GHz channel. Snapshots of the Dame et al. 
(2001) /=1^0 data are also shown in Fig. 7 but will be dis- 
cussed in Sect. 8.1. 

Dust emission is quite intense in both Orion and Taurus 
(Fig 7, right column) and, as can be seen in Fig. 6, some cor- 
relation is indeed found between the Planck CO maps and the 
dust template in these locations not unexpectedly. This is par- 
ticularly true in Orion where the images show very similar pat- 
tern of CO and dust. In the Polaris region the correlation loosens 
a lot; dust is far less intense and shows a diff'erent distribution 
than CO. For the weakest dust emission, probably in atomic gas, 
there is no CO signal. These correlation plots against dust should 
be compared with the ones performed against ground-based CO 
measurements, where the correlation is this time much tighter 
(this is discussed in Sect. 8.1, see Fig. 10). Looking at several 
other molecular clouds (e.g. Ophiucus, Chameleon), we find that 
whatever the method, the correlation with CO is always much 
tighter than the one with dust, making us confident that dust is 
not a major issue for CO-rich regions. 

Nonetheless, type 2 and type 3 maps do suff'er some level 
of dust contamination. This can be seen by eye in the CO(l^O) 
images of the Taurus regions in Fig. 7 (top) where the type 2 map 
shows more flux than the type 1 CO at the locations where the 
dust is the brightest. However, it is in the Galactic plane of the 
C0(2^1) TYPE 2 map that contamination is the most important. 
This is particularly visible in Fig. 4. 

6.2.2. Other CO lines 

For all of the Planck reconstructed CO maps we expect con- 
tamination from other CO isotopologue transition lines. Lines 
from the ^^CO isotopologue are the main contaminants in the 
TYPE 1 maps, while ^^CO and the ^^CO(3^2) both contribute to 



^ The temperature has little impact on the conversion coefficients of 
the dust. 
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Type 1 - Orion 
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Fig. 7. Images (12°. 5 x 12°. 5) of the Taurus (left column), Orion (middle column) and Polaris (right column) regions considered 
throughout the paper. The first row corresponds to the Dame et al. (2001) survey. The two middle rows show the type 1 and type 2 
/=1^0 Planck CO products. The Planck-Wl 545 GHz channel map, used as a dust template, is shown in the last row. All maps 
are smoothed to 30'. 
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the TYPE 2 maps. The contribution of each contaminant line can 
be estimated to first order as followed: 



TYPE 1 : Forgetting about other contaminants, the CO content of 
the TYPE 1 Planck map can be written as 



CO: 



typel 



'CO 



A^bolos A^bolos 



(5) 



where X represents the /=1^0, /=2^1 or J=3^2 transition, 
Wi are the weights of the linear combination and F\^^^ is the 

^^'^^CO transmission in bolometer /. Using the weights of the 
TYPE 1 linear combination and using the bandpass coefficients 
for ^^CO, we find the quantity 2f'°'" to be equal to 0.53 

at 115 GHz, 0.01 at 230 GHz and 0.36 at 353 GHz. Assuming 
a ^^CO to ^^CO ratio of 0.2, this translates in COl'^^^ ^ 1.11 x 

^^CO^-o and CO^^^^ ^ 1.07 x ^^CO^-^ i.e. an overestimation of 
about 10% of the type 1 /=1^0 and /=3^2 maps. The eff'ect 
on the /=2^1 transition is negligible. 



TYPE 2: The type 2 maps are constructed as C0^pg2 = Zv 
where X stands for /=1^0 or /=2^1, My is the Planck map at 
the channel frequency v and Wy is a weighting factor. Using the 
average CO transmissions in each channels, the CO content of 
the TYPE 2 maps formally reads 



co;., = "CO . ^^^^ 



^2CO''|l+W353<Fi2)353 



{Fn) 



(6) 



where (/^i2,i3) are the average conversion coeflScients of 
^^CO and '^^CO at 100 or 217 GHz, depending on the 
transition under scrutiny and <i^i2)353 is that of ^^CO 
in the 353 GHz channel. In the following, we assume 
the the line ratios to be ^^C0{?>^2) I^^C0(\^{)) = {)2 and 
^^CO{3^2)I^^CO{2^\) = QA as obtained from the median of 
the line ratio distributions computed from the type 1 maps, (see 
Planck Collaboration PIP 2013). The widths of the distributions 
are large and thus we consider uncertainties of ±0. 1 in the two 
ratios. Notice that within these uncertainties the values con- 
sidered here are consistent with those derived from the Ruler 
analysis in Sect. 4.2.3. Using these values we find C0{|~^^ ^ 



1.02 xi^co^^-'^ and COl^y^^ 



0.75 X i^CO^^ ^\ From this, we 
see that contamination by the /=3^2 line can be neglected in 
the /=1^0 TYPE 2 map, but reduces by ~ 25% the ^^CO sig- 
nal in the /=2^1 map. For this reason, we decided to correct 
the TYPE 2 map from this effect by dividing the output type 2 
{J =2 I) map by 0.75. With this correction, the overall calibra- 
tion of the TYPE 1 and type 2 (/=2 ^1) maps are in agreement 
(see Sect. 7). 

The second term of the contamination due to ^^CO can now 
be estimated. Using the bandpass estimate of the coefficients, the 
quantity {Fuqq)I{Fi2qq) is equal to 1.14 at 115 GHz and 0.82 at 
230 GHz, which translate into COf^p^^^^ ^ 1.2 x ^^CO^^-^^ and 
^^type2 ~ 1-2 X ^^CO^^"^^ when assuming a ratio of 0.2 between 
^^CO and ^^CO and the 0.75 correction factor for the /=2^1 
line. 



type2 



The values above are given as a rough estimate of the ^^CO 
contamination. The latter depends on the isotopic ratio and 
relative optical depth, and will, in practice, vary across the sky. 



6.2.3. Masking SZ clusters and point sources 

Planck CO maps are also contaminated by other localized non- 
CMB signals that have not been included in the component 
separation. This is the case of radio point sources, which are 
present in all three types of maps. Also present in the maps is 
the Sunyaev-Zeldovich (SZ) effect is a secondary anisotropy of 
the CMB coming from the interaction of CMB photons with the 
hot electron population of galaxy clusters. The Planck mission 
has shown the potential of galaxy cluster detection via the SZ 
effect with the publication of a catalogue of a few thousands 
clusters (Planck Collaboration XXIX 2013). As far as CO ex- 
traction is concerned, SZ is yet another contaminant foreground 
of the TYPE 2 and type 3 maps, as these rely on multi-channel 
information. On the other hand type 1 maps are not affected. 

For this first release of the Planck CO maps, a point source 
mask is provided with the type 1 and a point source -i- SZ mask 
for the TYPE 2 and type 3maps. The point source mask corre- 
sponds to Planck-YiF\ official 100 GHz mask (we chose the 4- 
cr, S/N = 10 mask) where we unmasked any pixels located 
within \b\ < 1.5°. This modification was necessary given that 
many point sources in the mask are molecular cold cores lo- 
cated in the Galactic disk, so that using the original mask hides 
most of the Galactic CO. For the SZ mask we start from a sim- 
ulated Compton parameter map of the Planck cluster sample 
(Planck Collaboration XXIX 2013; Planck Collaboration XXI 
2013) where a universal generalized pressure profile (Arnaud 
et al. 2010) was assumed. The mask is then generated by im- 
posing a threshold of 4 x 10"^ in Compton parameter units to 
this map. 



6.3. Absolute calibration uncertainties 

Three main contributors to the absolute calibration uncertainties 
of the Planck CO maps have been identified: i) the Planck-HFl 
calibration uncertainties of the temperature maps; ii) the uncer- 
tainty on the CO bandpass conversion coeflftcients; and iii) the 
uncertainty on the ^^CO contribution. We do not include the ef- 
fect of dust in this absolute calibration given that, conversely to 
^^CO, dust is not always spatially correlated to CO and is there- 
fore considered as a systematic effect (see Sect. 6.2.1). 

As described in ?, the absolute calibration of the Planck 
maps is about 0.55% from 100 GHz to 217 GHz, 1.25% at 
353 GHz, and 10% at 545 and 857 GHz. Uncertainties on the av- 
eraged bandpass CO conversion coefficients are at the 5% level 
at 100 GHz and 1% at 217 and 353 GHz (Planck Collaboration 
IX 2013). We described how ^^CO affected the Planck CO maps 
in Sect. 6.2.2; using the results obtained with ^^CO/^^CO = 0.2 
(Solomon et al. 1979), and assuming a conservative uncertainty 
of ±0.1 on this ratio, we can estimate the ^^CO contribution to 
the absolute uncertainty in each CO map. 

Final calibration uncertainties are summarized in Table 4. 
A conservative estimate of the total absolute calibration uncer- 
tainties of each map is given in the last column of the table. 
For most maps, the calibration uncertainties are of the order of 
10 % and are dominated by the ^^CO contribution. Notice that 
for TYPE 1 C0(2^1) and CO(3^2) the uncertainties are signifi- 
cantly smaller as they are less contaminated by ^^CO. 
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Table 4. Absolute calibration uncertainties of Planck CO maps. 



Type 




Line 


HFI calib. 


<^co> 


i^CO 


Total 








[%] 


[%] 


[%] 


[%] 


TYPE 


1 


CO(l^O) 


0.55 


5 


5 


10 


TYPE 


1 


C0(2^1) 


0.55 


1 




2 


TYPE 


1 


CO(3^2) 


1.25 


1 


3 


5 


TYPE 


2 


CO(l^O) 


0.70 


5 


10 


15 


TYPE 


2 


C0(2^1) 


0.70 


1 


10 


11 


TYPE 


3 




3 


5 
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7. Internal validation of the CO maps 

Before validating these maps on external data, it is possible to 
perform some internal checks by comparing the CO maps of dif- 
ferent types. As mentioned previously, the type 1 and type 2 CO 
products can differ in the Galactic plane because of the higher 
level of contamination affecting the type 2 products. However, 
at high Galactic latitudes, where no significant free-free or dust 
emission is expected, the products may be compared in order 
to assess the overall inter-calibration of the type 1, 2 and 3 CO 
maps. To perform this comparison, the maps are first smoothed 
to a common resolution of 30' and degraded to A^side = 128 to 
avoid noise correlation between samples. After masking point 
sources, correlation plots are produced for all remaining pixels 
located at Galactic latitudes \b\ > 25°. 

Figure 8 compares both type 1-type 2 and type 1-type 3 
maps. At high latitudes, CO emission is very sparse so that most 
pixels have very low emission. Therefore, below an empirical 
threshold, shown as the dashed line in each figure, the number 
of points is such that a contour representation of binning suffices. 
Above the dashed line, individual pixels are plotted directly and 
used to compute the best fit. The orange line gives the best-fit 
values and is compatible within errors with y = x. This indicates 
that the different types of product share the same overall calibra- 
tion, be it for the /=1^0 maps (Fig. 8, left and middle panel) or 
/=2^1 maps (Fig. 8, right panel). 

The different types of product may also be compared in spe- 
cific molecular clouds, as shown in the top and bottom panels of 
Fig. 9 for the /=1^0 (top) and /=2^1 (bottom) lines, respec- 
tively. As for the comparison to dust, the Taurus (left), Orion 
(middle) and Polaris (right) molecular clouds have been chosen 
because of the three very different molecular environments they 
host. Figure 9 shows tight correlations between the type 1 and 
TYPE 2 CO(l^O) products in these three molecular clouds. The 
best fits are computed for all points in the range [2-15] K km s"\ 
where the bulk of the emission lies. The fit uncertainties are 
dominated by the errors bars (not shown in the figure) of the 
TYPE 1 points. The type 2 CO(l^O) map shows more flux than 
the TYPE 1 with best-fit slopes of about 1.1 for Taurus, Orion and 
Polaris. Such behavior is nonetheless expected, given that the 
TYPE 2 CO(l^O) map suffers from more ^-^CO and dust contami- 
nation than the type 1 CO(l^O) map as described in Sects. 6.2.2 
and 6.2.1. Similar trends are observed for the C0(2^1) maps 
(bottom panel), but here we also give as dashed lines the best- 
fits that would have been obtained for a change of ±0.1 in the 
^^CO(3^2)/ ^^C0(2^1) ratio assumed to correct the type 2 
map from the ^^CO(3^2) contamination (see Sect. 6.2.2). Given 
the absolute calibration and statistical uncertainties quoted above 
for each of the maps, the residuals between the type 1 and type 2 
maps are consistent within the errors. 




5 10 

TYPE 1 CO [1-0] (Kkms-^) 




2 4 6 

type 1 CO [2-1] (Kkms-^) 



Fig. 8. Correlation between the type 1 and type 2 (top) and 
TYPE 1 and TYPE 3 (middle) /=1^0 CO maps at high Galactic 
latitudes 25° < \b\ < 90°. Correlation between the type 1 and 
TYPE 2 /=2^1 CO maps is shown on the bottom plot.The solid 
line corresponds to j = x. Below the dashed horizontal line the 
log colour scale gives the number of pixels in a given bin. The 
best linear fits (solid red line) are computed from the individ- 
ual pixels plotted above the dashed line. The best-fit slopes and 
uncertainties are: (a) (0.93 ± 0.17); (b) (0.87 ± 0.20); and (c) 
(0.96 ± 0.12). The intercept is compatible with zero within the 
error bars. 
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Taurus ^ Taurus 




TYPE 1 CO [1 -0] (K km s TYPE 1 CO [2-1 ] (K km s 



Fig. 9. Correlation between type 1 and type 2 CO maps at 1 15 GHz (left column) and 230 GHz (right column) in the Taurus (top), 
Orion (middle) and Polaris (bottom) molecular clouds. The best linear fits (orange solid line) have been computed between 2 and 
15 Kkm s"^ at 1 15 GHz and between 1 and 15 Kkm s"^ at 230 GHz, where bulk of the data lie and avoiding the noise around zero. 
The best-fit slopes and uncertainties are: (a) (1.12 ± 0.03); (b) (1.14 + 0.01); (c) (1.10 ± 0.04); (d) (1.07 + 0.01); (e) (1.07 ± 0.06); 
and (f) (1.17 ± 0.04). The intercepts are compatible with zero within error bars. The orange dot-dashed lines in the bottom panels 
correspond to the recalibration uncertainty due to the contribution of the J =3^2 line in the type 2 /=2^1 product. See text for 
detail. 



8. Comparison with external data 



In this section, we validate the Planck CO products described 
in Sect. 5 using existing ground-based CO data of the first three 
transition lines. 



8.1. Comparison with the Dame et ai. (2001) CO J =1^0 
survey 

The publicly available data of Dame et al. (2001) represent the 
most complete survey of Galactic ^^CO /=1^0 emission to 
date. The original data consist of a composite map constructed 



14 



Planck collaboration: Galactic CO emission in Planck 



from a set of 37 independent surveys taken by the 1.2 meter 
Millimetre- Wave Telescope at the CfA^. In order to compare 
the Planck CO maps to the Dame et al. (2001) data, we use 
the velocity-integrated HEALPix rendition of the survey which 
is available on the Lambda website^. 

8.1.1. Molecular clouds 

We focus again on the three molecular clouds Taurus, Orion and 
Polaris. Images of the Dame et al. (2001) and Planck data in 
these fields are shown in Fig. 7. The maps were smoothed to a 
common resolution of 30' and degraded to /"/side = 128 (to avoid 
noise correlation between samples) before plotting the correla- 
tions shown in Fig. 10. In each panel of this figure, the three 
types of Planck CO [1-0] maps (shown in diff'erent colours) are 
correlated with the Dame et al. (2001) data. The slopes of the 
best-fit linear regressions for the three products are given in the 
legend of each panel (they are found with a l-cr statistical error 
smaller than 1%). The spread in the data points is reminiscent of 
the level of noise in the maps, showing once again that the type 1 
map is the noisiest while the type 3 map has the best signal-to- 
noise ratio. 

The TYPE 1 CO is in good agreement with the Dame et al. 
(2001) data, with a slope of about 1.05 in both Taurus and Orion. 
The excess with respect to the 1-to-l correlation is at the level 
expected from ^^CO contamination (see Sect. 6.2.2). For Polaris 
the slope is ~ 1.1, but in such faint regions, the level of noise 
of the TYPE 1 map makes the comparison more difficult. For all 
clouds, the type 2 CO shows roughly a 20% excess compared to 
the Dame et al. (2001) data that is due to a combination of both 
dust and ^^CO contamination. It performs particularly poorly in 
Orion where the dust is significantly correlated with CO emis- 
sion. Finally, the type 3 CO does as well as the type 1 product in 
Taurus, but also suff'ers from some contamination in Orion with 
a correlation of 1.16. This test was conducted in other molecu- 
lar regions (not shown here) and similar trends were found. The 
level of the residuals is always compatible with the calibration 
uncertainty of the maps. 

From these, we conclude that the type Imap is the most ro- 
bust /=1^0 map in terms of CO extraction, but its low S/N 
ratio makes it unsuitable for the study of faint CO regions. The 
TYPE 2 and type 3 maps both suff'er from some level of dust con- 
tamination but allow us to probe fainter regions and to perform 
discovery studies. 

8.1 .2. Velocity effect in the type 1 /=1 ^0 map 

After checking the behaviour of our maps in these specific loca- 
tions, we focus here on the Galactic plane as a whole. Still work- 
ing at 30', the Planck map has been degraded to A^side = 512 
in order to match the HEALPix resolution of the Dame et al. 
(2001). The Planck map has been recalibrated to the Dame et al. 
data by performing a linear regression on all pixels (COtypei ~ 
1 . 16x ^^CODame), after which operation the diff'erence was taken. 
Figure 11 shows the residual between the type 1 CO(l^O) map 
and Dame et al. (2001) Galactic plane composite survey. We ob- 
serve a clear modulation pattern in the Galactic disk, produc- 
ing a positive residual at longitude Igai < 90° and a negative 
residuals at Igai > 270°. This residual can reach 5% of the to- 
tal intensity of CO and it is interpreted as the signature of the 

^ The complete data set can be retrieved from http://www.cfa. 
harvard . edu/rtdc/CO/ 
^ http://lambda.gsfc.nasa.gov/ 
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Dame et al, CO [1-0] (Kkm s"^) 

Fig. 10. Correlation between the three types of CO products with 
the Dame et al. (2001) data in the Taurus (top), Orion (middle) 
and Polaris (bottom) molecular clouds. The best-fit slope and 
intercept are given in Table 5. Notice that the intercepts are in all 
cases below 0.7 K km s"^ . 



rotation of the Galactic disk. Because the CO line frequency is 
shifted due to the Doppler eff'ect, the emission along diff'erent 
lines-of-sight with diff'erent velocities produce a diff'erent CO re- 
sponse in a given bolometer. The eff'ect is then averaged over all 
bolometers in the final reconstructed CO map. Using the mean 
velocity of CO emission in the Dame et al. (2001) survey, we 
can fit the residuals as a linear function of the velocity. We find 
that this eff'ect at the /=1^0 transition is well represented by 
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Table 5. Best linear fit parameters, slope and intercept, for the correlation plots shown in Fig. 10 



Panel type 1 type 2 type 3 



(a) (1.054 + 0.016); (-0.231 + 0.097) (1.195 + 0.001); ( 0.409 + 0.070) (1.081 + 0.001); ( 0.334 + 0.003) 

(b) (1.056 + 0.009); ( 0.141 + 0.069) (1.290 + 0.001); (-0.654 + 0.005) (1.163 + 0.001); (-0.488 + 0.002) 

(c) (1.104 + 0.029); (-0.070 + 0.052) (1.156 + 0.002); ( 0.102 + 0.004) (1.012 + 0.001); ( 0.112 + 0.002) 




Fig. 11. Residual map with respect to the Dame et al. (2001) map 
for the TYPE 1 /=1^0 Planck CO map. The Planck maps have 
been recalibrated to Dame before performing the subtraction. 



^^typel = + ^/^OO ^"0 ^^fme^^' ^^^^^ ^ ^^^^^^ 

locity of the gas. This formula was also found consistent with the 
diff'erence observed between survey maps (for which the Planck 
satellite velocity is modulated by the satellite revolution around 
the Sun) and may be used to correct the Planck map for this ve- 
locity eff'ect. Note that this eff'ect is not seen when performing 
the residuals in the same way for the type 2 and type 3 producta 
as it is i) somewhat averaged out by the use of channel maps in- 
stead of bolometer maps; and ii) hidden by the higher Galactic 
plane contamination of these maps. 

8.2. Comparison with NANTEN II ^^CO and ^^COJ=1^0 
data 

Observations in both ^^CO and ^^CO have been carried out with 
the NANTEN II millimetre-submillimetre telescope at Atacama, 
between February 2010 and October 2012. The 4-m dish pro- 
vides a half-power beam width (HPBW) of 2.'6 in ^^CO /=1^0. 
The spectrometer is a digital Fourier spectrometer with a fre- 
quency resolution of 61 kHz and a bandwidth of 1 GHz. A 
Hamming window function is applied when the FFT is per- 
formed in the spectrometer, which results in a final frequency 
resolution of 79.3 kHz. The velocity coverage and resolution are 
~ 2600 kms-i and 0.21 kms'^ for ^^CO, and 2720 kms"^ and 
0.22 kms"^ for ^^CO respectively. All observations have been 
carried out in On The Fly (OTF) mode, where the telescope con- 
stantly drives across the aree, with an output grid of 60''. The 
standard size of an OTF block is 1° by 1° and at least one scan 
of longitudinal direction and lateral direction have been done 
for each OTF block. Each scanning data set was combined by 
the basket weaving method (Emerson & Graeve 1988) to re- 
duce scanning noise in each OTF direction. The intensity cal- 
ibration was made with the chopper-wheel method (Kutner & 
Ulich 1981) and absolute intensity calibration performed on p- 
Oph, Ori-KL, M17SW, and Perseus. 

The soon-to-be-published data we use here consist of a 9° x 
2° Galactic plane patch around Galactic coordinates (315°, 0°). 
We smoothed the NANTEN and Planck maps to a 15' resolution 
before carrying out the comparison. This resolution is chosen 
as it corresponds to that of the type 2 CO map and allows a 
significant noise reduction in the type 1 map (which comes with 
a 9. 65 'native resolution). Remembering that we expect a 0.53 x 



^^CO contribution to the type 1 map (see Sect. 6.2.2), we show 
in Fig. 12 the combined NANTEN II ^^CO -h 0.53 x ^^CO on the 
top panel, the type 1 map in the middle and the correlation plot 
in the bottom panel. The agreement between the two sets of data 
is good with a correlation of 0.95. While not shown here, we 
perform this test for the type 2 CO map as well, using this time 
1.14 X ^^CO as is estimated in Sect. 6.2.2. The agreement is once 
again satisfactory with a best-fit slope of 1.00. In both cases, the 
level of the residuals between the best-fit and the combination of 
NANTEN data is found to be < 10% and thus compatible with 
the calibration uncertainties. 

This test is important, as it validates our understanding of the 
^^CO contamination in the Planck CO(l^O) maps. It also high- 
lights that the type 2 CO(l^O) map does not suff'er from ma- 
jor dust contamination in the Galactic plane since, had this been 
the case, the correlation with NANTEN would have presented a 
slope greater than one. We also performed an independent check 
of the ^^CO contamination based on FCRAO data, which corrob- 
orates the results obtained with NANTEN. This is the purpose of 
the next section. 

8.3. Comparison with FCRAO ^^CO J=1^0 data 

Along the same lines as the comparison performed with 
NANTEN II data, it is possible to gain further confidence 
in our estimation of the ^^CO contamination using FCRAO 
^^CO(l^O) data along the ^2CO(1^0) data of Dame et al. 
(2001). In this analysis, we use the publicly available Boston 
University FCRAO Galactic Ring Survey"^ (FCRAO-GRS) that 
has been described in Jackson et al. (2006). It consists of a 
75.4 deg^ survey of the Galactic plane, between 18°< / <55?7 
and \b\< r at 46" resolution. 

We first re-project the FCRAO-GRS data into the HEALPix 
pixelization scheme (A^side = 2048) using a nearest grid-point ap- 
proach and correct for the 0.48 beam efficiency given in Jackson 
et al. (2006). Binning over the ratio ^^COfcrao/^^COd ame? WC 
compute the average ratio <COtypei)/(^^CODame) in each bin and 
plot the correlation in Fig. 13. Error bars are obtained assuming 
white noise, the amplitude of which is estimated at high Galactic 
latitude (\b\ > 60°). The standard deviation in each bin is then 
simply estimated using the number of pixels in the bin. We as- 
sume the noise in the Dame et al. (2001) map to be negligible. 

The correlation between the two ratios is good and from the 
best-fit of the linear regression we see that the Planck type 1 map 
can be written as 

CO[yp"'i^ = 0.53 X ^^COfcRAO + 1.07 X ^^COoame , 

which is agreement with the 0.53 x^^ CO contamination estimate 
we made from the bandpass coefficients (see Sect. 6.2.2). This 
test shows, in a completely independent manner to that of the 



9 The data of the BU-FCRAO-GRS may be retrieved at http : // 
www . bu . edu/Galacticring/ 
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Fig. 12. Comparison of the Planck (/ = 1 ^ 0) type 1 map to 
the NANTEN II survey in a 9° x 2° section of the Galactic plane. 
Top: Combined NANTENII map as ^^CO + 0.53 x ^^CO, where 
the factor applied to ^^CO corresponds to the bandpass estimate 
of the ^^CO contamination in the type 1 map (see Sect. 6.2.2). 
Middle: Planck type 1 /=1^0 map. Bottom: correlation plot be- 
tween the two maps. The log colour scale represents the number 
of pixels in a given intensity bin. The solid and dashed lines 
represent y = x and the best linear fit, respectively. The best- 
fit slope and intercept are (1.007 + 0.001) and (5.243 ± 0.032), 
respectively. 



NANTEN comparison, that we have an excellent grasp on the 
CO content of the type 1 map. 

8.4. Comparison with the AMAN0GAWA-2SB CO J=2^ 1 
survey 

The AMAN0GAWA-2SB survey (Handa et al. 2012; Yoda et al. 
2010) carried out simultaneous Galactic plane (/ = 10° - 245°) 
observations of the ^^CO and ^^C0(2^1) transitions with the 
Tokyo-NRO 60-cm telescope, with a resolution and grid-spacing 
of 9 and 3.75' respectively. Here, we use a 35° x 8° section of 
the AMAN0GAWA-2SB ^^CO Galactic plane survey for com- 
parison to the Planck C0(2^1) maps. 

The area under scrutiny is shown in Fig. 14 for the 
AMANOGAWA data set (top) and the Planck type 1 CO [2-1] 
map (middle) at the working resolution of 15'. The correlation 
between the two data sets in given in the bottom panel and shows 
remarkable agreement with a best-fit value for the correlation of 
0.95. Conversely to the CO(l^O) map, no ^^CO contribution 
was required to match the two data sets. This is once again in 



Best fit 




tf.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 

^^COpcRAO / ^^COoame 

Fig. 13. Average ratio <COtypei)/(^^CODame) in bins of the 
^^COpcRAo/^^CODame ratio. The best linear fit is represented 
by the orange solid line. The best-fit slope and intercept are 
(0.53 ± 0.03) and (-1.07 ± 0.01), respectively. 



agreement to what was predicted in Sect. 6.2.2 where we found 
that the ^^CO contribution to the type 1 /=1^0 line should be 
negligible, given the weighted bandpass transmission of the iso- 
topologue in the 217 GHz channel. The residuals between the 
best-fit and the AMANOGAWA data are found to be about 5%, 
which is compatible with the combined calibration uncertainties 
of Planck and AMANOGAWA. 

8.5. Comparison with the FIR AS CO J =1^0, J =2^1 and 
J =3^2 surveys 

Fixsen et al. (1999) have analyzed a by-product of the 
Far Infrared Absolute Spectrophotometer (FIRAS on COBE 
launched in 1989) sky maps, namely integrated velocity maps 
of well-known Galactic lines. These include carbon, nitrogen 
and water lines, but, more importantly for the present study, 
the complete ^^CO rotation ladder. These are absolutely cali- 
brated all- sky measurements. FIRAS has a better spectral res- 
olution (13.5 GHz) than Planck-HFl but a much lower angu- 
lar resolution. FIRAS integrated CO lines are retrieved from the 
Lambda website^^. The first three transitions at 115.27, 230.54 
and 345.80 GHz come from the low frequency emission line 
maps. They are converted from nWm"^sr"^ to Kkms"^ with 
the coefficient 10"^^c^/2/:y^Q. Error maps are processed in the 
same way and we neglect gain uncertainties. 

Planck CO maps are convolved with the FIRAS beam in the 
same way as described in Appendix A of ? and compared to 
FIRAS line maps within the FIRAS pixelization scheme. This is 
simply done by taking the Planck value at the grid point nearest 
to the FIRAS pixel within 15'. Uncertainties are much larger for 
FIRAS than Planck; the noisiest FIRAS values are not kept in 
the comparison. The cut-off' is chosen as 160, 1 1 and 7 Kkm s"^ 
for the CO(l^O), C0(2^1), and CO(3^2) lines, respectively. 
The FIRAS Galactic centre values may not be reliable because 
of the large velocity spread and we therefore discard any pixels 
within 4° of the Galactic centre. There is no possible tuning in 
this straightforward comparison. 



http : //lambda . gs£c . nasa . gov 
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Fig. 14. Comparison of the Planck J=2-l type 1 CO map to the 
AMAN0GAWA-2SB Galactic plane survey. From top to bottom 
we display the AMANOGAWA data, the type Idata in the same 
region and the correlation between the two maps. The log colour 
scale represents the number of pixels in a given intensity bin. 
The solid orange and dashed blue lines represent y = x and the 
best linar fit, respectively. The best-fit slope and intercept are 
(0.941 ± 0.001) and (0.792 ± 0.009), respectively. 



Fig. 15 shows the correlation plot between the type 1 Planck 
CO data and the FIRAS /=1^0 (top), /=2^1 (middle) and 
/=3^2 (bottom) data. In each panel, the dashed blue line cor- 
responds to J = X, while the orange line gives the best fit. We 
find that FIRAS and Planck all-sky CO data are consistent with 
a linear relationship. For the CO(l^O) line, Fixsen et al. (1999) 
claim a detection only in the Galactic centre. With Planck we ob- 
tain a statistical detection of CO(l^O) with FIRAS outside the 
Galactic centre at the 7 cr level. For the two other lines, FIRAS 
measurements tend to overestimate the CO emission with re- 
spect to Planck by at most 10%. This is well within the absolute 
calibration error for Planck and beam uncertainties for FIRAS. 



8.6. Comparison with HARP/ACSIS CO J =3^2 data 

Surveys in of the CO(3^2) transition are not as numerous as 
those of lower transitions, which makes assessing the quality of 
the TYPE 1 /=3^2 map more difficult. We used part of the JCMT 
HARP/ACSIS i2CO(3^2) Galactic plane data (soon to be pub- 
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Fig. 15. Comparison of the type 1 Planck CO maps with the 
FIRAS /=1^0 (top), /=2^1 (middle) and /=3^2 (bottom) 
data. The solid orange and dashed blue lines represent y = x and 
the best linear fit, respectively. The best-fit slopes and intercepts 
are: (a) (0.98 ± 0.17) and (0.33 ± 0.65); (b) (1.12 ± 0.02) and 
(-0.01 + 0.04); and (c) (1.13 ± 0.03) and (0.04 ± 0.02). 



lished). Details on the HARP/ACSIS system and calibration can 
be found in Buckle et al. (2009); Polychroni et al. (2012). 

We use a 2° long section of the Galactic plane from / = 12.5 
to / = 14.5 where the survey has latitude range ±0.5°. The orig- 
inal HARP data were not clipped at some n cr threshold but col- 
lapsed over the whole spectral range, after performing a heavy 
smooth on all regions with signal < 3cr. This ensures there is 
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Fig. 16. Correlation plot between the latitude-collapsed longitu- 
dinal profiles (see text for details) of type 1 /=3^2 HFI CO map 
and the HARP J =3 ^2 survey. The solid orange and dashed blue 
lines represent y = x. . 



no significant emission missing, while minimizing the addition 
of noise. The data were then smoothed to about the 4.5' resolu- 
tion, to mimic the Planck resolution but with pixels of around 
5''. The two maps were then brought down to the same pixeliza- 
tion, the least constraining rebinning resulting in a 43 x 23 pixels 
patch, with about 2.7' pixels. Because of the noise level inherent 
to the TYPE 1 map, a large dispersion is seen when simply per- 
forming the correlation between all pixels of these two patches. 
Therefore, we perform another rebinning, collapsing the latitude 
dimension to get a 43 point "average" longitudinal profile of 
the two maps. The correlation of these two profiles is plotted 
in Fig. 16. 

Despite our rebinning, the correlation remains quite noisy, 
but the figure shows that the y = x line is a fair representation 
of the data when ignoring the two outliers. More data would be 
required to perform a thorough characterisation of the type 1 
/=3^2 map. 

8. 7. High Galactic latitude structures 

In order to achieve a deep, high-Galactic-latitude statistical com- 
parison of the Planck CO maps with ground-based observa- 
tions, we use the observations performed at CfA and described 
in Hartmann et al. (1998) (HMT98) and Magnani et al. (2000) 
(MHHSTOO) for the North and South Galactic hemisphere re- 
spectively. These data consist of 8' beam observations of the 
\b\ > 30° sky, on a regularly-spaced grid, observable from 
Cambridge (MA). In total, 10,443 points in the north Galactic 
hemisphere and 4,934 in the south were observed. A detection 
threshold of 0.3 KRj-km-s"^ (3cr) was considered. In total 26 de- 
tections were reported in the North and 133 in the South. 

First, for each of these detection points, we compute the cor- 
responding Planck TYPE 2 and type 3 CO [1-0] fluxes (in a 10' 
diameter disc). The type 1 map was not selected for this analysis 
as its noise level makes it unsuitable for high Galactic latitudes. 
The top panel in Fig. 17 shows the correlation of the Planck 
TYPE 2 and type 3 fluxes as a function of the 159 HMT98 and 
MHHSTOO fluxes. A reasonable agreement and compatible ab- 
solute calibration is reached in these faint CO-emitting regions. 
The bottom panel presents the histogram of all three types of 
Planck CO maps fluxes at the 15,218 locations where no detec- 




-2 2 

Planck CO [1-0] flux (K km s'^) 



Fig. 17. Comparison of the Planck CO(l^O) maps with the 
high Galactic latitudes CO detections and non detections from 
MHHSTOO. Top: correlation between the CO detections in 
MHHSTOO and the Planck type 2 and type 3 maps. The lat- 
ter was first smoothed to 15'. Bottom: distribution of all Planck 
CO(l^O) maps (at 15' resolution) at the non-detection positions 
of MHHSTOO. The width of the distribution is compatible with 
the noise in the maps. Diff'erences with respect to the values in 
Table 3 can be explained by the inhomogeneities of the noise in 
the Planck maps. 



tion is expected from the CfA surveys. The Gaussian distribu- 
tions obtained are all centred on zero and have a standard devi- 
ation compatible with the noise level of the map. Discrepancies 
with respect to the values in Table 3 can be explained the inho- 
mogeneities of the noise in the Planck maps. This comparison 
to the HMT98 and MHHSTOO highlights the high quality of the 
Planck CO maps in the faintest CO regions. 



9. Conclusions 

The Planck-HFl maps contain a non-negligible component that 
is attributed to the emission, in the HFI spectral bands, of the 
three first CO rotation transitions coming from the interstellar 
medium. Several component separation methods are used to iso- 
late this component from the others (CMB and diff'use Galactic 
emission, mainly thermal dust), by using the spectral diversity 
and large coverage of Planck-HFl. 
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This paper presents all-sky maps of the velocity-integrated 
emission of CO in the /=1^0, /=2^1, and J=3^2 transitions. 
Three different types of maps are produced: i) type 1 maps ob- 
tained from a single-channel analysis for the three transitions; ii) 
TYPE 2 maps available for the first two transitions only and com- 
ing from a multi-channel analysis; and iii) a single type 3 map 
relying on a multi-line approach. 

We have characterized the maps in terms of resolution, noise and 
systematic eff'ects by a cross-comparison and internal validation 
of the maps. Comparison to external CO data set was also per- 
formed and a general agreement is found, always within the es- 
timated calibration uncertainties of the data. Note however that 
the TYPE 1 J=3^2 map could not be as thoroughly tested as the 
others because of the few external CO(3^2) data sets available. 

In summary, type 1 maps are the noisiest but suff'er from 
the least systematic errors and foreground contamination. The 
TYPE 2 CO products provide high S/N single-line maps at the 
price of larger systematics and foreground contamination, and 
are of lower resolution. Finally, the type 3 map gives a high S/N 
and high resolution CO map but relies on the assumption of con- 
stant line ratios across the sky. Both type 2 and type 3 maps are 
more suitable for a use at intermediate and high-Galactic lati- 
tudes. The TYPE 1 maps are, by contrast, more reliable in terms 
of foreground contamination in the Galactic plane and on high 
CO intensity regions. In all these maps the ^^CO contribution is 
not corrected and needs to be accounted for in any further scien- 
tific analysis. 

The Planck CO maps are used within the Planck collabora- 
tion for the definition of regions of the sky where CMB, thermal 
SZ and cosmic infrared background cosmological studies can 
safely be performed. The Planck CO maps presented in this pa- 
per are delivered as part of the Planck nominal mission public 
archive. 

As discussed in Planck Collaboration PIP (2013), Planck 
multi-line CO survey reveals that the bulk of the Galactic CO 
emission is subthermally excited and originates in a diff'use com- 
ponent. The highest sensitivity CO maps prove that CO emission 
extends well beyond the known boundaries of molecular cloud. 
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Appendix A: Description of the Planck CO products 

The characteristics of the released maps are the following. We 
provide HEALPix pixelization all-sky maps with Nside=2048. 
For one transition, the CO velocity-integrated line signal map is 
given in K.km units at the CO transition frequency. A conver- 
sion factor from this unit to the native unit of HFI maps (Kcmb) 
is provided in the header. Four maps are given in a single file 
per transition (1-0), (2-1) and (3-2) and per type 1, type 2, and 

TYPE 3, 

- a CO signal map; 
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- a CO null test map; 

- a standard deviation map; 

- a mask for which a value of 1 corresponds to a valid pixel. 

TYPE 1 products have the native HFI resolution 9 .'65, 4.'99 
and 4.'82 for the CO [1^0], [2^1] and [3^2] transitions re- 
spectively. TYPE 2 products have 15' resolution and are given for 
the two first transitions. The type 3 product has a 5. '5 effective 
resolution and corresponds to a combined map of the three CO 
transitions normalized at the (/ = 1 ^ 0) transition frequency. 
These maps are shown in the left column of Figs., 3, 4 and 5. 
The corresponding standard deviation maps are displayed in the 
right column. 

Appendix B: Validation of CO extraction methods 
on simulations 

Planck component separation methods (see Planck 
Collaboration XII 2013) are generally validated using the 
FFP6 set of simulations (Planck Collaboration ES 2013), which 
aim at providing a complete realization of the Planck mission. 
They rely of the Planck Sky Model (PSM) fully described in 
?. FFP6 maps can be used for the multi-channel and multi-line 
approaches, but they do not provide the individual bolometer 
maps required to validate the single-channel CO extraction 
method. For the latter, CO-tailored simulations are needed, as 
described below. 



We apply the multi-channel CO extraction method described 
in Sect. 4.2.2 to the set of FFP6 intensity maps and characterize 
the reconstructed CO(l^O) map with the input map of the FFP6 
simulations. As for some of the other tests done throughout the 
paper, the comparison is performed in the Orion, Taurus, and 
Polaris molecular clouds. 



A linear correlation (C0type2 = ^ COffpg + b) between the 
two maps is found in the three regions, with slopes of 0.99+0.01, 
1.013 ± 0.001, 0.999 ± 0.002 for Taurus, Orion, and Polaris re- 
spectively. We also checked the standard deviation of the resid- 
ual maps C0ffp6 - C0type2 against the average cr in each re- 
gion obtained from the type 2 CO(l^O) standard deviation map 
(see Fig. 3). In Polaris, the residual is within the ±\-cr level, 
while it is at about 1.5-(T for Taurus and Orion. In these two 
latter regions, the slightly poorer behaviour of the residual is at- 
tributed to the simple modified blackbody dust modelling of the 
multi-channel approach that cannot properly capture the dust of 
the FPP6 simulations. With these caveats, the multi-channel ap- 
proach is successfully validated on FFP6 simulations. 



B.1. Single-channel approach: TYPE 1 maps 

We use here specific simulations of the Planck bolometer maps 
at 100 GHz. For each bolometer in the Planck 100 GHz channel, 
we construct a map of the sky emission considering CMB and 
foreground components (Galactic diff'use emission from syn- 
chrotron, thermal dust and free-free, and Galactic and extra- 
Galactic point sources). We use the same Planck Sky Model 
(PSM) model as for the FFP6 simulations (see below for a de- 
tailed description). For the CO emission we just use the Dame 
et al. (2001) map as a template and applied the CO bandpass con- 
version coeflftcients from Table 2 to convert to Kcmb units. The 
per bolometer null maps were added to account for the noise 
contribution. The MILCA algorithm was then applied to these 
maps to reconstruct the CO signal using the single channel ap- 
proach discussed in Sect. 4.2.1. We found that when the ^^CO 
conversion coefficients are perfectly known the CO emission is 
well reconstructed. The contribution from foregrounds was esti- 
mated to be less than 0.5 % of the CO signal overall and the rms 
of the residuals consistent with the noise at about the l-cr level. 

We also tested the sensitivity of the MILCA algorithm to un- 
certainties on the CO conversion coefficients. These uncertain- 
ties translate into noisier final CO maps and an overall calibra- 
tion bias. For 1% and 5 % uncertainties on the CO conversion 
coefficients we observe a combined signal to-noise-reduction of 
6% and 60%, respectively. 

B.2. Multi-channel approach: type 2 maps 

The multi-channel approach may be validated using Planck 
FFP6 set channel intensity maps. Relevant to our purpose here, 
it includes in particular a y"^-^^ free-free emission component, 
a two-component dust emission from model 7 of Finkbeiner 
et al. (1999), and CO emission using the CO map of Dame et al. 
(2001) as the template. 



B.3. Multi-line approach: type 3 map 



The quality of the reconstruction of the type 3 CO map was 
also validated using the FFP6 simulations. As described in 
Sect. 4.2.3, the type 3 CO map is a byproduct of the CMB 
oriented component separation procedure and therefore, valida- 
tion on simulations is presented in the companion CMB -oriented 
component separation paper Planck Collaboration XII (2013). 



APC, AstroParticule et Cosmologie, Universite Paris Diderot, 
CNRS/IN2P3, CEA/lrfu, Observatoire de Paris, Sorbonne Paris 
Cite, 10, rue Alice Domon et Leonie Duquet, 75205 Paris Cedex 
13, France 

^ Aalto University Metsahovi Radio Observatory, Metsahovintie 114, 
FlN-02540 Kylmala, Finland 

^ African Institute for Mathematical Sciences, 6-8 Melrose Road, 
Muizenberg, Cape Town, South Africa 

^ Agenzia Spaziale Italiana Science Data Center, c/o ESRIN, via 
Galileo Galilei, Frascati, Italy 

^ Agenzia Spaziale Italiana, Viale Liegi 26, Roma, Italy 

^ Astrophysics Group, Cavendish Laboratory, University of 
Cambridge, J J Thomson Avenue, Cambridge CB3 OHE, U.K. 

^ Astrophysics Research Institute, Liverpool John Moores University, 
Twelve Quays House, Egerton Wharf, Birkenhead CH41 ILD, U.K. 

^ Astrophysics & Cosmology Research Unit, School of Mathematics, 
Statistics & Computer Science, University of KwaZulu-Natal, 



21 



Planck collaboration: Galactic CO emission in Planck 



Westville Campus, Private Bag X54001, Durban 4000, South Africa 



California, Los Angeles, CA 90089, U.S.A. 



^ CITA, University of Toronto, 60 St. George St., Toronto, ON M5S 
3H8, Canada 

10 CNRS, IRAP, 9 Av. colonel Roche, BP 44346, F-31028 Toulouse 
cedex 4, France 

11 California Institute of Technology, Pasadena, California, U.S.A. 

1^ Centre for Theoretical Cosmology, DAMTP, University of 
Cambridge, Wilberforce Road, Cambridge CB3 OWA U.K. 

1^ Centro de Estudios de Fisica del Cosmos de Aragon (CEFCA), 
Plaza San Juan, 1, planta 2, E-44001, Teruel, Spain 

Computational Cosmology Center, Lawrence Berkeley National 
Laboratory, Berkeley, California, U.S.A. 

1^ Consejo Superior de Investigaciones Cientificas (CSIC), Madrid, 
Spain 

16 DSM/Irfu/SPP, CEA-Saclay, F-91 191 Gif-sur-Yvette Cedex, France 

1^ DTU Space, National Space Institute, Technical University of 
Denmark, Elektrovej 327, DK-2800 Kgs. Lyngby, Denmark 

1^ Departement de Physique Theorique, Universite de Geneve, 24, 
Quai E. Ansermet,1211 Geneve 4, Switzerland 

1^ Departamento de Fisica Fundamental, Facultad de Ciencias, 
Universidad de Salamanca, 37008 Salamanca, Spain 

^0 Departamento de Fisica, Universidad de Oviedo, Avda. Calvo 
Sotelo s/n, Oviedo, Spain 

^1 Department of Astronomy and Astrophysics, University of Toronto, 
50 Saint George Street, Toronto, Ontario, Canada 

Department of Astrophysics/IMAPP, Radboud University 
Nijmegen, RO. Box 9010, 6500 GL Nijmegen, The Netherlands 

Department of Electrical Engineering and Computer Sciences, 
University of California, Berkeley, California, U.S.A. 

Department of Physics & Astronomy, University of British 
Columbia, 6224 Agricultural Road, Vancouver, British Columbia, 
Canada 

Department of Physics and Astronomy, Dana and David Dornsife 
College of Letter, Arts and Sciences, University of Southern 



^6 Department of Physics and Astronomy, University College London, 
London WCIE 6BT, U.K. 

Department of Physics, Gustaf Hallstromin katu 2a, University of 
Helsinki, Helsinki, Finland 

Department of Physics, Nagoya University, Chikusa-ku, Nagoya, 
464-8602, Japan 

Department of Physics, Princeton University, Princeton, New 
Jersey, U.S.A. 

^0 Department of Physics, University of California, One Shields 
Avenue, Davis, California, U.S.A. 

^1 Department of Physics, University of California, Santa Barbara, 
California, U.S.A. 

Department of Physics, University of Ilhnois at Urbana-Champaign, 
1110 West Green Street, Urbana, Illinois, U.S.A. 

Dipartimento di Fisica e Astronomia G. Galilei, Universita degli 
Studi di Padova, via Marzolo 8, 35131 Padova, Italy 

Dipartimento di Fisica e Scienze della Terra, Universita di Ferrara, 
Via Saragat 1, 44122 Ferrara, Italy 

Dipartimento di Fisica, Universita La Sapienza, P. le A. Moro 2, 
Roma, Italy 

^6 Dipartimento di Fisica, Universita degli Studi di Milano, Via 
Celoria, 16, Milano, Italy 

Dipartimento di Fisica, Universita degli Studi di Trieste, via A. 
Valerio 2, Trieste, Italy 

Dipartimento di Fisica, Universita di Roma Tor Vergata, Via della 
Ricerca Scientifica, 1, Roma, Italy 

Discovery Center, Niels Bohr Institute, Blegdamsvej 17, 
Copenhagen, Denmark 

^'^ Dpto. Astrofisica, Universidad de La Laguna (ULL), E- 3 8206 La 
Laguna, Tenerife, Spain 

"^1 European Space Agency, ESAC, Planck Science Office, Camino 
bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo, 
Villanueva de la Canada, Madrid, Spain 

European Space Agency, ESTEC, Keplerlaan 1, 2201 AZ 
Noordwijk, The Netherlands 

1^ Finnish Centre for Astronomy with ESO (FINCA), University of 
Turku, Vaisalantie 20, FIN-21500, Piikkio, Finland 

^ Graduate School of Science and Engineering, Kagoshima 
University, 1-21-35 Korimoto, Kagoshima, Kagoshima 890-0065, 



22 



Planck collaboration: Galactic CO emission in Planck 



Japan 



Grenoble, F-38041, France 



Helsinki Institute of Physics, Gustaf Hallstromin katu 2, University 
of Helsinki, Helsinki, Finland 

INAF - Osservatorio Astrofisico di Catania, Via S. Sofia 78, 
Catania, Italy 

INAF - Osservatorio Astronomico di Padova, Vicolo 
deir Osservatorio 5, Padova, Italy 

INAF - Osservatorio Astronomico di Roma, via di Frascati 33, 
Monte Porzio Catone, Italy 

^'^ INAF - Osservatorio Astronomico di Trieste, Via G.B. Tiepolo 11, 
Trieste, Italy 

5^ INAF/IASF Bologna, Via Gobetti 101, Bologna, Italy 

^1 INAF/IASF Milano, Via E. Bassini 15, Milano, Italy 

INFN, Sezione di Bologna, Via Imerio 46, 1-40126, Bologna, Italy 

INFN, Sezione di Roma 1, Universita di Roma Sapienza, Piazzale 
Aldo Moro 2, 00185, Roma, Italy 

IPAG: Institut de Planetologie et d'Astrophysique de Grenoble, 
Universite Joseph Fourier, Grenoble 1 / CNRS-INSU, UMR 5274, 



lUCAA, Post Bag 4, Ganeshkhind, Pune University Campus, Pune 
411 007, India 

Imperial College London, Astrophysics group, Blackett Laboratory, 
Prince Consort Road, London, SW7 2AZ, U.K. 

Infrared Processing and Analysis Center, California Institute of 
Technology, Pasadena, CA 91125, U.S.A. 

Institut Neel, CNRS, Universite Joseph Fourier Grenoble I, 25 me 
des Martyrs, Grenoble, France 

Institut Universitaire de France, 103, bd Saint-Michel, 75005, Paris, 
France 

Institut d'Astrophysique Spatiale, CNRS (UMR8617) Universite 
Paris-Sud 11, Batiment 121, Orsay, France 

61 Institut d'Astrophysique de Paris, CNRS (UMR7095), 98 bis 
Boulevard Arago, F-75014, Paris, France 

6^ Institute for Space Sciences, Bucharest-Magurale, Romania 

6^ Institute of Astronomy and Astrophysics, Academia Sinica, Taipei, 
Taiwan 

Institute of Astronomy, University of Cambridge, Madingley Road, 
Cambridge CB3 OHA, U.K. 

6^ Institute of Astronomy, University of Tokyo, 2-21-1 Osawa, Mitaka, 
Tokyo, Japan 

66 Institute of Theoretical Astrophysics, University of Oslo, Blindern, 
Oslo, Norway 

6^ Instituto de Astrofisica de Canarias, CA^ia Lactea s/n. La Laguna, 
Tenerife, Spain 

6^ Instituto de Fisica de Cantabria (CSIC-Universidad de Cantabria), 
Avda. de los Castros s/n, Santander, Spain 

6^ Jet Propulsion Laboratory, California Institute of Technology, 4800 
Oak Grove Drive, Pasadena, California, U.S.A. 

Jodrell Bank Centre for Astrophysics, Alan Turing Building, School 
of Physics and Astronomy, The University of Manchester, Oxford 
Road, Manchester, M13 9PL, U.K. 

^1 Joint Astronomy Centre, 660 N. Aohoku Place, University Park, 
Hilo, Hawaii 96720, U.S.A. 

Kavli Institute for Cosmology Cambridge, Madingley Road, 
Cambridge, CB3 OHA, U.K. 



73 



LAL, Universite Paris-Sud, CNRS/IN2P3, Orsay, France 



LERMA, CNRS, Observatoire de Paris, 61 Avenue de 
I'Observatoire, Paris, France 

Laboratoire AIM, IRFU/Service d'Astrophysique - CEA/DSM - 
CNRS - Universite Paris Diderot, Bat. 709, CEA-Saclay, F-91191 
Gif-sur-Yvette Cedex, France 

^6 Laboratoire Traitement et Communication de 1' Information, CNRS 
(UMR 5141) and Telecom ParisTech, 46 rue Barrault F-75634 Paris 



23 



Planck collaboration: Galactic CO emission in Planck 
Cedex 13, France U.S.A. 

Laboratoire de Physique Subatomique et de Cosmologie, Universite University of Granada, Departamento de Fisica Teorica y del 
Joseph Fourier Grenoble I, CNRS/IN2P3, Institut National Cosmos, Facultad de Ciencias, Granada, Spain 
Poly technique de Grenoble, 53 rue des Martyrs, 38026 Grenoble 

cedex, France University of Miami, Knight Physics Building, 1320 Campo Sano 

Dr., Coral Gables, Florida, U.S.A. 

Laboratoire de Physique Theorique, Universite Paris-Sud 11 & 

CNRS, Batiment 210, 91405 Orsay, France Warsaw University Observatory, Aleje Ujazdowskie 4, 00-478 



Warszawa, Poland 



Lawrence Berkeley National Laboratory, Berkeley, California, 
U.S.A. 

Max-Planck-Institut fur Astrophysik, Karl-Schwarzschild-Str. 1, 
85741 Garching, Germany 

McGill Physics, Ernest Rutherford Physics Building, McGill 
University, 3600 rue University, Montreal, QC, H3A 2T8, Canada 

MilliLab, VTT Technical Research Centre of Finland, Tietotie 3, 
Espoo, Finland 

National Astronomical Observatory of Japan, 2-21-1 Osawa, 
Mitaka, Tokyo 181-8588, Japan 

Niels Bohr Institute, Blegdamsvej 17, Copenhagen, Denmark 

Nobeyama Radio Observatory, National Astronomical Observatory 
of Japan, 462-2 Nobeyama, Minamimaki, Minamisaku, Nagano 
384-1305, Japan 

Observational Cosmology, Mail Stop 367-17, California Institute of 
Technology, Pasadena, CA, 91125, U.S.A. 

Optical Science Laboratory, University College London, Gower 
Street, London, U.K. 

SB-ITP-LPPC, EPFL, CH-1015, Lausanne, Switzerland 

SISSA, Astrophysics Sector, via Bonomea 265, 34136, Trieste, Italy 

School of Physics and Astronomy, Cardiff University, Queens 
Buildings, The Parade, Cardiff, CF24 3AA, U.K. 

Space Research Institute (IKI), Russian Academy of Sciences, 
Profsoyuznaya Str, 84/32, Moscow, 117997, Russia 

Space Sciences Laboratory, University of California, Berkeley, 
California, U.S.A. 

Special Astrophysical Observatory, Russian Academy of Sciences, 
Nizhnij Arkhyz, Zelenchukskiy region, Karachai-Cherkessian 
Republic, 369167, Russia 

'^^ Stanford University, Dept of Physics, Varian Physics Bldg, 382 Via 
Pueblo Mall, Stanford, California, U.S.A. 

Sub-Department of Astrophysics, University of Oxford, Keble 
Road, Oxford OXl 3RH, U.K. 



96 



Theory Division, PH-TH, CERN, CH-121 1, Geneva 23, Switzerland 



9^ UPMC Univ Paris 06, UMR7095, 98 bis Boulevard Arago, 
F-75014, Paris, France 

9^ Universite de Toulouse, UPS-OMP, IRAP, F-31028 Toulouse cedex 
4, France 

Universities Space Research Association, Stratospheric Observatory 
for Infrared Astronomy, MS 232-11, Moffett Field, CA 94035, 



24 



